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Abstract 

Classical Jacobi polynomials P^ 1 ^, with a,0 > — 1, have a number of well-known 
properties, in particular the location of their zeros in the open interval (—1,1). This 
property is no longer valid for other values of the parameters; in general, zeros are com- 
plex. In this paper we study the strong asymptotics of Jacobi polynomials where the real 
parameters a n , f3 n depend on n in such a way that 

lim — = A , lim — = B , 

n — >oo Ji n — >oo ji 

with A, B EM.. We restrict our attention to the case where the limits A, B are not both 
positive and take values outside of the triangle bounded by the straight lines ^4 = 0, 
B = and A + B + 2 = 0. As a corollary, we show that in the limit the zeros distribute 
along certain curves that constitute trajectories of a quadratic differential. 

The non-hcrmitian orthogonality relations for Jacobi polynomials with varying pa- 
rameters lie in the core of our approach; in the cases we consider, these relations hold 
on a single contour of the complex plane. The asymptotic analysis is performed using 
the Deift-Zhou steepest descent method based on the Riemann-Hilbert reformulation of 
Jacobi polynomials. 



Corresponding author. 
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1 Introduction 

Jacobi polynomials Pi a '^ are given explicitly by 

{z) = 2 -n ±( n n + _ a V n + P) (*-!)*(* + 1)- * , (1.1) 
fc=0 ^ ' ^ ' 

or, equivalently, by the Rodrigues formula [321 Ch. IV] 

P^(z) = ~ l)- a (z + 1)"" [(z - l) n+a (z + l)"+»] . (1.2) 

In the classical situation (a, /3 > —1) the Jacobi polynomials are orthogonal in [—1,1] with 
respect to the weight function (1 — x) a (l + x)^ and, consequently, their zeros are simple and 
located in (—1, 1). 

Expressions 1)1.1(1 - 1)1.21) show that the definition of P^ may be extended to arbitrary 
a,/3 G M (or even C); but some properties of the classical polynomials, in particular the 
location and simplicity of the zeros, are no longer valid. In fact, P^*'^ may have a multiple 
zero at z = 1 if a G {—1, • • • , — n}, at z = — 1 if (3 G {—1, • • • , — n} or at z = oo (which means 
a degree reduction) if n + a + (3 G {—1, • • • , —n}. 

More precisely, for k G {1, . . . , n}, we have (see [221 formula (4.22.2)]), 

p(-fc,fl(^_ nn + P+l) (n-k)\ (z-l\ k (ktP) 



This implies in particular that P„ (z) = if additionally maxjfc, — 0} < n < k — (3 — 1. 

Analogous relation; 
k + I < n, we have 



Analogous relations hold for P^ a ' ^ when / G {l,...,n}. Thus, when both t,I 6 N and 



Pt k '- l) (z) = 2~ k ~ l (z - l) k (z + 1)' ffl, (*) • (1-4) 
Furthermore, when n + a + /3 = —A: G { — 1, • • • , — n}, 

see jSH Eq. (4.22.3)]; we refer the reader to [321 §4.22] for a more detailed discussion. Taking 
into account formulas l)1.3|) - ()1.5(l we exclude these special integer parameters from our further 
analysis. 

In this paper we study the asymptotic behavior of Jacobi polynomials P^ n ^ n \ where the 
parameters a n , f3 n depend on the degree n in such a way that 

lim^L=A, lim^ = 5, (1.6) 

n— too 11 n— >oo n 
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and 

A<0<B, A + B>-1, A^-l, (1.7) 

In [21], the authors considered different regions of the (A, U)-plane, corresponding to different 
cases in the asymptotic study of Jacobi polynomials with varying parameters satisfying (|1.6|) . 
The symmetry relations (see |32) Ch. IV]) 

P^\z) = (-l)"ff '<*)(*) (1.8) 

and 

P(^)(z) = (i^) n p(^)(i±|), (1.9) 

where a' = —2n — a — (3 — 1, allow to restrict our study to the following cases, from which 
all the others can be obtained: 

A,B>0; (C.l) 

A<-1, and A + B>-1; (C.2) 

-1<A<0 and B>0; (C.3) 

A + B>-1, and A,B<0; (C.4) 

A + B<-1 and A,B>-1 (C.5) 

(see Figure ^ which appeared first in |24j . where equivalent regions under those transforma- 
tions are shown). 

Case C.l is classical and has been widely studied (see [2], [3], [HI [22] and [2E])- 

The asymptotic results therein are based on either the well known orthogonality conditions 
satisfied by the Jacobi polynomials on [—1, 1], or on their integral representation. 

However, until very recently, strong asymptotics of p^ a "'^ n ^ when a n , (3 n take arbitrary 
real values and limits 1)1.6)) exist, has not been established. In this case the orthogonality 
conditions were unknown and the complex saddle points make the application of the classi- 
cal steepest descent method to the integral representation of Jacobi polynomials practically 
unfeasible. 

A non-hermitian orthogonality satisfied by Jacobi polynomials in the case C.2 has been 
observed in [21]; this fact was used there to establish the asymptotic zero distribution using 
a potential theory approach. 

Recently in 23 a whole spectrum of orthogonality conditions for Jacobi polynomials 
with arbitrary real parameters has been established. In particular, we can find examples of 
orthogonality on a contour or arc of the complex plane, an incomplete or quasi-orthogonality, 
or even multiple or Hermite-Pade orthogonality conditions. The classification of the cases 
depends on the number of inequalities — 1 < A < 0, — 1 < B < 0, —2 < A + B < — 1 that are 
satisfied. In particular, cases C.3, C.4 and C.5 correspond to combination of parameters A 
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A=-1 A=0 



B=0 



B=-1 



A+B=-1 



A+B=-2 



Figure 1: The five different cases in the classification of Jacobi polynomials with varying 
parameters according to |24"j . 

and B such that exactly one, exactly two, or exactly three, respectively, of the inequalities 
are satisfied (cf. Figure P). 

Nevertheless, the method used in |24j cannot be immediately extended to the rest of 
the cases. One of the essential assumptions there is a non-hermitian orthogonality of the 
polynomials on a single contour, on which the support of a certain equilibrium measure has 
a connected complement. 

Due to this reason, in the steepest descent method of Deift and Zhou (see based 
on a matrix Riemann-Hilbert problem, was used to establish the strong uniform asymptotics 
of the Jacobi polynomials with parameters satisfying conditions C.5. In this paper, we use 
several results and ideas from there. 

The aim of the present article is to extend this analysis to sequences of Jacobi polynomials 
with varying parameters corresponding to cases C.2-C.3. Note that along with the case C.l, 
these are the only situations when a full system of orthogonality relations on a single contour 
in C exists. 

We also remark that a similar study, but for Laguerre polynomials with varying parame- 
ters, has been carried out in |19|-|20| and |25j . 

The paper is organized as follows. In Section 2, the main results about strong and 
weak zero asymptotics are formulated, along with some preliminary definitions and lemmas 
which are proved in Section 3. In Section 4, a full set of orthogonality relations on a single 
contour allows to pose a Riemann-Hilbert problem and to apply the Deift and Zhou's steepest 
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descent method (see (3]) to transform the Riemann-Hilbert problem in order to obtain strong 
asymptotics of the solution. Finally, the last section is devoted to the proofs of the main 
results. 



2 Main results 
2.1 Basic definitions 

Let us denote by C+ = {z £ C : lm(z) > 0} and C" = {z G C : lm(z) < 0}. For A,B 
satisfying (|1.7j) . we define the points 



B 2 - A 2 ± 4V(A + l)(B + 1)(A +^TT) 
Cl ' 2 = (A + B + 2) 2 • (2J) 

We will use the following convention: for (A,B) such that A < — 1 < A + B (case C.2), 
Ci G C+ and C2 = Ci; for -1 < A < < B (case C.3), we agree that -1 < Ci < C2 < 1- 
With these £1,2 > consider the set 

N *Z jz £ C : Re jT {{t - dt = o| , (2.2) 

where we continue the integrand analytically along the path of integration. Obviously, the 
set does not depend on the branch of the square root. In fact, it coincides with the union of 
the critical trajectories of the quadratic differential 



(*-Cl)(*-C2) ,2 



{Z 2 - l) 2 



dz z , (2.3) 



or more precisely, with their projection on C. Taking into account the local structure of 
trajectories of quadratic differentials (see e.g. [2B] or [HH)> we can prove the following (see 
Figure HI): 

Lemma 1 If parameters A, B satisfy condition (|1.7|) . then for defined in (12.1(1 the 

quadratic differential is regular. In other words, all its critical trajectories are finite and have 
the following global structure (see Figure^) : 

• For (A, B) such that A < —1 < A + B (case C.2), N consists of three arcs which 
connect ("1,2 and intersect the real line in exactly one point, in such a way that each of 
the intervals (— 00,— 1), (—1,1), (1, 00) is cut by one of these arcs. 

• For (A,B) such that —1<A<0<B (case C.3), N consists of three arcs; one of 
them is the real interval [Ci 5 C2] and the other two are Jordan contours, passing through 
z = (1 (respect., z = Qi) and enclosing z = — 1 (respect., z = 1). 
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Figure 2: Typical structure of the set N for cases C.2 (left) and C.3. Dashed lines on the 
left are orthogonal critical trajectories. 

Now we define some relevant curves. We denote by V the rightmost curve from N . For 
case C.2, V consists of an arc connecting £12 and crossing once the interval (1, +00), and for 
case C.3, it is a closed contour passing through z = (2 and surrounding z = 1. For case C.2 
we also consider the orthogonal trajectories AA^ (defined by replacing Im in (|2.2|) by Re). 
As in Lemma n if is easy to prove that their global structure is as appears in Figure |2J left 
(dashed lines). We denote by 7 + the arc of Af ± joining £1 and —1, and 7" =7+. 

Finally, we define the set E as the smallest connected subset of M containing £12 and T. 
Namely, 



As we see, in case C.2 the set E is made of one critical trajectory of the quadratic differential 
()2.3j) . while in case C.3 it is made of two. In both cases E is oriented from £j to £21 and, in 
case C.3, clockwise, in such a way that (1, +00) is cut from the upper to the lower half-plane 
(see Figure |3J). For any function / analytic and single- valued in a neighborhood of E, this 
selection of the orientation induces two boundary values of / on E that we denote by /+ and 
/_ depending if we approach E from the left or from the right, respectively. On the sequel, 
we shall make use of the concept of the polynomial convex hull of E, which is denoted by 
Pc(E). In case C.2, Pc(E) = T, so that Int(Pc(E)) = 0, where by Int(e) we denote the set of 
inner points of e. Analogously, in case C.3, Pc(E) is the union of S and of the closure of the 
bounded component of its complement, given by Int(Pc(E)). 

Next we define some functions that will play a role in what follows. We denote 




if A < -1< A + B, (case C.2) ; 
if -1<A<0<B, (case C.3) . 



(2.4) 



R(z) d = f ^{z-Ci)(z-C2). 
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-1 • 




Figure 3: Contours E for cases C.2 (left) and C.3. 



It is a multi- valued and analytic function in C, and we select its single- valued branch in a 
plane cut from £x to £2 by imposing that 

z^oo ^ 

This allows us to define the (a priori complex) measure: 

o/x(2j = — _ 2 tte, (2.5) 

with E defined in (|2.4j) and oriented as explained. By Lemma [2] below, fj, is a unit positive 
measure on E. Using Cauchy's Theorem it is easy to find an analytic expression for its 
Cauchy transform: 

« rdrtt) = A+B+2 m _An_Bj2 L , 6CWE); (2 . 6) 

PV 1 J z-t 2 z 2 -lz-lz + i' \v;> v; 

additionally, in case C.3, 

, , A + B + 2 R(z) A/2 B/2 _ .„ „. 

= 2^^H-^T-ITT' ^Int(Pc(E)). (2.7) 

Now, let us a define in C \ E a function which plays a key role in the description of the 
strong asymptotics of Jacobi polynomials, 

G(z) = exp ^ dtj . (2.8) 
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We normalize G by imposing that 



lim G{z) = 1 , (2.9) 



where the limit is taken with z approaching £2 from C\T (in case C.2) or from C + \E (in case 
C.3). Observe that since jl is the Cauchy transform of a unit measure on E, function G is 
analytic and single- valued in C \ E in both cases considered. Taking into account l|2.6|) — (|2.7|) . 
we see that there exists 



«^Iim^. (2.10) 



In addition, let 



W (z)=w(z;A,B)^c(z-l) A /\z + l) B ^=e^^ + <*«) , (2.11) 

which is a multi- valued analytic function in C\{±1}. In what follows, we fix its single- valued 
analytic branch in C\ (— 00, 1], by choosing the constant c (or the path of integration) in such 
a way that 

lim 10(2) = 1 , (2.12) 

where again the limit is taken with z approaching (2 from C\T (in case C.2) or from C + \ E 
(in case C.3). 

The last ingredient for the asymptotics is given by the functions 

Nll {z)^ a{z)+a{zyl and N 12 (z) = a{z) ~ a{Zrl (2.13) 
2 2% 

(this notation is chosen because they will be entries of a certain matrix N, see (14.15(1 ). where 

1 

a{z) ={^zT (2 - i4) 

is defined inC\rfor^<-l<A + J B (case C.2), and in C \ [&, C2] for -1 < A < < B 
(case C.3). We select the branch of a imposing the normalization condition 

lim a(z) = 1 . 
Then, Nn(z) — > 1 and N\2(z) -^Oasz^ 00. 
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2.2 Strong asymptotics 

First, we consider the strong asymptotics for Jacobi polynomials with varying parameters 
satisfying Q1.6)l - (|1.7|) with z away from E. 

Theorem 1 Let (A,B) satisfy (|1,7|) . Then, for n — * oo, the monic Jacobi polynomials 
Pn = p( An ' Bn ^ have the following asymptotic behavior: 

Pn(z)=(^y N n (z) (l + o(~\\ , (2.15) 

locally uniformly in C \ Pc(£), where constant k was defined in (|2.1()jl . 
Furthermore, in the bounded component ofC \ X, 



Pn (z) =JL ((G(^ 2 (^)) n ]\ru(*) (i + o Q 

+2ie -An7ri sin ( A7rn ) G n (z)iYi 2 (^) (l + O 



I 

n 



(2.16) 



In particular, this theorem shows that zeros of p^ An ' Bn ^ do no t accumulate in C \ Pc(S). 
Next, we describe the asymptotics on S, but away from the branch points (j,2 : 

Theorem 2 Let (A,B) satisfy ([1.7)1 . Then, for n —* oo, i/ie monic Jacobi polynomials 
p n = p( An ' Bn ^ have the following asymptotic behavior for z away from £x 2 ■ 



(a) In case on £/te "±"-side ofT 

(2.17) 



ft>(*) ~ (g"(z)JV u (z) (l + oQ 



± (G(^ 2 (z)piV 12 (z) (i + oQ)) . 



(b) In case C.5, on the "—"-side ofT formula (|2.16|) is still valid, while on the "+"-side of 

r, 

Pn (z) =-L (G n {z)N ll (z)(l + 0^- 



-2ie- Anwi sin 



/ / x xxx ( 2 - 18 ) 

(Avrn) (G(z)u; 2 (z)) n N 12 (z) (l + O (- 
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(c) In case C.3, on the "+"-side of (Ci ? C2) 

Pn(z) ~ (&{z)Nu(z)(l + 0(l 



+ {G(z)w 2 (z)) n N 12 (z)(l + o(^ 



(2.19) 



while on the "—"-side 0/(^1,(2), 

Pn (z) =-L (G n (z)N u (z)(l + O f 1 



g— 2Amn 



n , , 

(2.20) 



{G{z)w 2 (z)) n Nl2 {z)( K l + o[^ 



Remark 1 All these asymptotic expressions match on the boundaries of the overlapping 
domains and on the respective regions. For instance, it will be shown that in a small neigh- 
borhood of every point of E (distinct from £1,2)1 \G{z)w{z)\ > 1 for z £ E. Hence, the first 
term in l|2.17j) - ()2.20j) is dominant, and away from E they reduce to ([2.15)1 . 

Furthermore, in the case C.3, if An are not exponentially close to integers (in the sense 
that will be made more precise below), the second term in p. 16)1 dominates, and we may 
write 

p n (z) = 2ie~ An ™ S m(Airn)N 12 (z) (l + O Q 

The asymptotic formulas above are no longer valid close to the branch points Ci,2- As 
it usually happens in a neighborhood of the "soft ends" of the support of an equilibrium 
measure, asymptotics is described in terms of the Airy function Ai(z) and its derivative. We 
give explicit formulas only for the rightmost (according to the orientation of E) branch point 
£2, which is in a certain sense, the "interesting" one. Clearly, the analysis at the other point 
is similar. 

In order to formulate our result in a more compact form it is convenient to introduce the 
function 

, , , def A + B + 2 f z R(t) 



dt. (2.21) 



Consider case C.2. Given a sufficiently small e > 0, and a neig hborhood A e (C 2 ) = {z € C : 
\ z ~ C2I < £ }i it is a single-valued analytic function in A e (^2) \ T. Furthermore, taking into 
account the local behavior of R, function 

/(z) d = f ^ (0(z)f 3 (2.22) 

can be extended as single- valued to the whole A £ (^2)- Here the 2/3rd power is chosen such 
that f(z) > along 7". 
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Theorem 3 Let (A,B) such that A < -1 < A + B (case C.2). Then, there exists e > 
such that if \z — C2I < we have that the monic Jacobi polynomials p n = p( An > Bn ^ satisfy: 

K n w n {z) \ a(z) \ \n 



a(z) 



^AiVVW) (i + o(i) 



where a(z) is defined in (|2.14|) . and we tofce f l / i {z) > along j~. 

Remark 2 Obviously, the asymptotic behavior near £1 in this case is completely symmetric 
to £2 with respect to R. 

Consider case C.3. For a sufficiently small < e < 1 — C21 <f> is single-valued and analytic 
in A e (C2) \ (Ci) C2)- Function /, defined again by formula (|2.22|) . can be extended as a single- 
valued function to the whole A e (^2), with the 2/3rd power chosen such that f(z) > along 

(C2,l). 

Theorem 4 Let (A,B) such that — 1 < A < < B (case C.3). Then, there exists e > 
such that if \z — C2I < £ , we have that the monic Jacobi polynomials p n = p^ An > Bn ) satisfy: 



' n l/6/l/4(s) 



K, n w n (z) \ a(z) \ \n J J 

-^k)^" 2 ' 3 ^ ( 1+ °©)) 



(2.23) 



where a{z) is defined in (|2.14|) . 

.4(t) = 4(4; A, n) = e"^ 4 ™ Ai(t) + 2ie™ /3 sin(^Trn) Ai (e 4 ™/ 3 t 
and we take / 1//4 (z) > along ((2, 1)- 

Remark 3 Formulas stated in Theorems are locally uniformly continuous both on the 
z and (A, B) planes, which allows to extend them to the general case of {a n } and {f3 n } 
satisfying (fTTS)) - (fl~7|) . 

2.3 Weak zero asymptotics 

As a corollary of the asymptotic formulas stated in the previous section we can obtain the 
distribution of the zeros of the sequence of polynomials P^* n, ^ n \ where {a n } and {/? n } satisfy 
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()1.6j) - ()1.7j) . By "weak zero asymptotics" we understand here the limit (in the weak-* sense) 
of the normalized zero counting measures associated with p^ n ^ n \ 

The measure [i introduced in (|2.5[) will be all we need for the description of the asymptotic 
behavior of the zeros in the case C.2. However, region C.3 comprises the pathological cases 
given by 1)1. 3|) . By continuity, we may expect here a variety of limit behaviors. In fact, in 
order to characterize completely the weak zero asymptotics of Jacobi polynomials in the case 
C.3 we need to use a 1-parametric family of measures including (|2.5|) , Namely, in case C.3, 
we must consider the sets 

M r =A/-j^) d = f {zeC: \G(z)w(z)\ = e r / 2 } =|.£C: RejT *gLdf = A + £ + 2 } ' 

for r > 0. They also consist of trajectories of the quadratic differential (|2.3|) . and Mq = M. 
Now, we define T r as the rightmost curve in J\f r or, what is the same, the part of J\f r which 
is entirely contained in the half-plane {z £ C : Hez > £2}- It is easy to check that for r > 
the level curve T r is a closed contour inside r = Tq surrounding the point z = 1 (see Figure 

n. 

0.4 r- 

0.3 - 

0.2 - 

0.1 - 
- 
-0.1 - 
-0.2 - 
-0.3 - 
-0.4 - 

-1 -0.5 0.5 1 

Figure 4: Some trajectories of the quadratic differential 1)2.3 J) , or equivalently, some level sets 
r r , for the values A = —0.8 and B = 0.5. 

For each r G [0, oo) we define the absolutely continuous measure 

= A+ 2 7Tl +2 T=^ dz ' *e(Ci,C2)ur rj (2.24) 
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and, for r = oo, the measure 



dfioo(z) = -A6i H — Xtft.Ca] « z ' ( 2 - 25 ) 

where X[Ci,Ca] * s ^ ne characteristic function of the interval [Ci ; C2] ; an d 5\ is the Dirac delta 
(unit mass point) at z = 1. 

Lemma 2 7/ (-4,5) satisfy (|1.7j) . i/ten /or r > measures \x r in (|2.24j) - ([2.25)1 ('and, in 
'particular, measure jj, in \2.5\) ) are unit positive measures. Moreover, for (A,B) satisfying 
-1 < A < < B (case C.3), we have that for < r < +00, 

M[Ci,C 2 ]) = l + A Mr(IY) = -A, (2.26) 
where we consider Tq = V and = {1}. 

Now we are ready to state the weak zero asymptotics for the Jacobi polynomials Pn an '^ n \ 
In the case C.3, when —1 < A < < B, we make an additional assumption: the sequence of 
parameters a n satisfies that the following limit exists: 

lim [dist(a n ,Z)] 1/n = e~ r , < r < +00. (2.27) 

n— >oo 

Then, it holds: 

Theorem 5 Consider a sequence of Jacobi polynomials P^ ln, ^ n \ n G N, such that sequences 
{a n }, {f3 n } satisfy M.b}) and \1.7\) . Then: 

(i) If (A,B) satisfy condition C.2, then the zeros of P r [ an ' /3n \ n G N, accumulate on the 
arc T, and measure n in \2. 5)) is the weak* limit of the corresponding normalized zero 
counting measures. 



(ii) If (A, B) satisfy condition C.3, and \2.27\) holds for some < r < +00, then the zeros 
of Jacobi polynomials Pn "' 13 "^, n G N ; accumulate on [Ci , C2] U T r; and measure fir- 
defined above is the weak* limit of the normalized zero counting measures. 

As we said before, part (i) of Theorem corresponding to case C.2, was established in 
|24| for parameters a n ,/3 n varying according to (jl.6f) but with A, B < —1, which is a region 
of the (^4, l?)-plane equivalent to A < — 1 < A + B by means of transformations (|1. 8(1 - 1)1. 9)1 . 

In the case C.3 the situation when r = 0, that is, 

lim [dist(a n , Z)] 1 ^™ = 1, 
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is generic, because it takes place when parameters do not approach the integers exponentially 
fast. When the limit in (|2.27|) is smaller than one, i.e. r > 0, curve T is replaced by a level 
curve T r , strictly contained inside the bounded component of the complement to T, and the 
support of the limit measure becomes disconnected. Finally, when parameters a n ,n S N, 
tend to integers faster than exponentially, the limit measure has a discrete part consisting of 
a Dirac mass at z = 1. We could have anticipated this phenomenon observing the coalescence 
of zeros given by (jl.3j) . 

Examples of zeros of Jacobi polynomials for Cases C.2 and C.3 are represented in Figure 

El 

Remark 4 The case C.3 deals with situations when real zeros arise. When — n < a n < 
and (3 n > -1, Pt ,P) satisfy a quasi-orthogonality relation (see Theorem 6.1 in [22]) which 
ensures the existence of, at least, n — [— a n ] zeros in (—1, 1). This lower bound of the number 
of zeros in (—1, 1) is exact, according to the so-called Hilbert-Klein formulas |321 Theorem 
6.72]. Since limn—nx) w ~[~ an l = 1 + A, looking at (|2.26|) we see that the mass of the part of 
the asymptotic measure of zeros supported on [Ci ; C2] C (—1, 1) agrees with the limit of the 
ratio of zeros placed in (—1, 1), given by the Hilbert-Klein formulas. 

Remark 5 At this point, it is natural to ask about what happens when A = —1 and B > 0, 
which is a transition case between C.2 and C.3. By (|1.<S|) - (|1.9|) . it describes also the situation 
when (^4, B) belongs to any of the straight lines A = — 1, B = — 1 and A + B = —1, outside 
of the square (^4, B) £ [—1, 0] x [—1, 0]. Roughly speaking, in this case the endpoints £1,2 are 
confluent in a single point, say £, and zeros approach a simple closed contour emanating from 
this point (critical trajectory) and surrounding z = 1, or other closed trajectories strictly 
contained in the interior of this critical trajectory. In [21], the equation of this critical 
trajectory is conjectured. This conjecture is proved in and ^2] for the particular case 
where a n = —n — 1 and j3 n = kn + 1, k being a fixed positive real number. This and the 
other transitions between cases C.l - C.5 deserve a separated treatment. 

3 Proof of the auxiliary lemmas 

Proof of Lemma ^ This proof is based upon the local structure of the trajectories of 
quadratic differentials (see [2S] or [SJ). We restrict our attention to the case where — 1 < 
A < < B (case C.3). The proof of the other case is similar (see also the proof of Lemma 
2.1 in [21])- 

First, we see that for — 1 < A < < B the quadratic differential (|2.3|) possesses two 
simple zeros at £12. Thus, we know that three critical trajectories emanate from £^2 at equal 
angles. Moreover, the segment [Ci , C2] C M, which is straightforward to verify by definition 
of M. 
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-1.5 -1 -0.5 



Figure 5: Up, case C.2: zeros of -Piqo 110 10 ' 5 ° 10 \ along with the curve T, correspond- 
ing to A = —1.1 and B = 0.5. Down, case C.3: zeros of P^ 80-10 ' 50 ~ 10 •* (left) and 
-P^ 80 10 ,5 ° 10 "* (right), along with the set S, corresponding to A = —0.8 and B = 0.5. 
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On the other hand, Q2.3j) has double poles at z = ±1 and z = do, in such a way that if we 
consider the rational function Q(z) = —(z — Ci)( z ~ (,2)(z 2 — l) -2 , we have that the residues 
of \/Q at these points are purely imaginary. Therefore, we conclude that near these double 
poles the trajectories are simple closed contours. 

Now, the symmetry of Q with respect to the real axis, along with the facts that the 
trajectories cannot tend to infinity and 1)2.3}) has no other singular point, allows one to ensure 
that the other critical trajectories are two closed contours emanating from (j and £2- The 
fact that a closed trajectory needs to surround a singular point implies that these closed 
trajectories intersect the real axis in two points, one of them in (l,+oo) and the other in 
(-oo,-l). 

□ 



Proof of Lemma 121 Taking into account the definition of R{z) and (|2.1j) . it is easy to see 
that: 

2A 

- < 0, if A < -1< A + B (case C.2) ; 

A + B + 2 ( 3-1 ) 

>0, i£-KA<0<B (case C.3) . 



R(l) 



2A 



{ A+B+2 



In the same way, 

R(-l) 



2B 
A + B + 2 



<0, if^<0<B and A + B>-1 



(3.2) 



Thus, the definition of £ in ()2.4I) yields that (|2.5|) is real- valued on £ and does not change sign 
on each of its components. The same remains valid when —1 < A < < B and < r < +00 
for n r on its support S r = T r U [Ci, C2]- 

Moreover, for (A, B) such that —1 < A < < B, taking into account the definition of E 
and S r , the residue theorem and (J3.1j) - (|3.2j) . we have for < r < +00: 



(J, r (T r ) 



J dnr(t) = -(A + B + 2) res (J^T[) = (A + -B + 2) 



R(l) 



clearly, also /ioo(roo) = /Uoo({l}) = -A. 
On the other hand, for < r < +00, 



Mr([Cl,C2]) 



C2 

dfj, r (t) 
Ci 

A + B + 2 
2 

A + B + 2 



res 



R(z) 



1 V z l - 1 
A 



+ res 



1 + 



i\z 2 -l 
B 

A+B+2 ~ A+B+2 



+ res 

z=oo 



Z 2 - 1 



1 + A, 
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and, therefore, 



l_i r (E r ) = / dfi r (t) = / dfi r (f) + / dfx r (t) = 1. 
Analogously, for (A, B) such that A < — 1 < A + B, it is easy to see that 



M (r) = / d M (t) = i , 




□ 



4 Riemann-Hilbert analysis 

4.1 Orthogonality and the Riemann-Hilbert problem 

As it was mentioned in the Introduction, the key fact for the asymptotic analysis is a full 
system of orthogonality relations satisfied by the Jacobi polynomials on simple contours, 
which allows to pose a matrix Riemann-Hilbert problem (RHP) and apply the Deift-Zhou 
steepest descent method. 

The following result was established in |231 Theorem 5.1]: 

Lemma 3 Let C be a Jordan arc connecting z = — 1 + Oi with z = — 1 — Oi and surrounding 
z = 1 once. If (3 > 0, then we have: 



where w(-;a,(3) has been introduced in (|2.11[) . 

From the seminal work of Fokas, Its and Kitaev ^1 (see also 0]) it is known that the 
orthogonality (|4.1|) can be characterized in terms of the following Riemann-Hilbert problem: 
find a matrix valued function Y : C\C — > C 2x2 satisfying the conditions below: 

(RH1.1) Y is analytic in C \ C. 

(RH1.2) Y has continuous boundary values on C, denoted by Y + and Y_, such that 




(4.1) 




for z G C. 



(RH1.3) As z -> oo, 
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(RH1.4) Y is bounded in a neighborhood of z = —1. 

Proposition 1 (|23j) The unique solution of the Riemann-Hilbert problem (RHl.l)-(RHl. 4) 

is given by 



Y(z) 



Qn-l(z) r^-Jc 



Pn (t)w 2 (t;a,f3) 
t-z 

gn -i(t)w 2 (t;q,/3) 



alt 
dt, 



where p n (z) = Pn (z) is the monic Jacobi polynomial, and q n -i(z) 
some suitable non-zero constant b n —^. 



h 1 P 



(«,/3), 



z), for 



Let (A,B) be a pair satisfying (|1.7|) . Then for every n £ N, monic polynomials p( An ' Bn } 
satisfy the conditions of Lemma |3J Hence, for any C as described in Lemma |3 polynomials 



P 
1 n 



(An,Bn) 



can be characterized as the (1, 1) entry of the unique matrix Y solving the Riemann- 



Hilbert problem (RH1.1)-(RH1.4) with 



a 



An, (3 = Bn. 



Taking advantage of the freedom in the selection of the contour C in ()4.1|) we will choose 
different configurations for both cases C.2 and C.3. This choice is mainly suggested by the 
numerical evidence on the actual location of zeros. For describing the appropriate C we will 
use the contours defined in subsection 12,11 corresponding to A and B fixed. 

In the case C.2, when A < —1 < A + B, we will make the contour C in Lemma 
coincide with 7~urU7 + , oriented clockwise (Figure|Hl left). Hence, we are interested in the 
Riemann-Hilbert problem (RH1.1)-(RH1.4) with C = 7" U V U 7+, a = An and j3 = Bn. 




Q. 




Figure 6: Contours C for the Riemann-Hilbert problem for Y; cases C.2 (left) and C.3. 



However, in case C.3, when —1 < A < < B, it is convenient to make part of the original 
contour C coalesce along the interval [—1, £2] (traversed twice in opposite directions), and the 
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rest go along the arc T. This deformation creates a new contour, which we denote again by 
C, and we choose its orientation as in Figure |HJ right. Consequently, it yields a new RHP, 
still characterizing the polynomials p n . With respect to the problem (RH1.1)-(RH1.4), we 
have to modify only the jump matrix on (— 1, C2) as a result of the coalescence of two original 



sub-arcs of C: its (1,2) entry becomes w 2 ^ ■ 
problem is: find a matrix valued function Y 
the following conditions hold: 

(RH2.1) Y is analytic on C \ ([-1, C2] U T). 



w 



2n 



on (—1,(2)- The new Riemann-Hilbert 
Y(A,B) : c \ ([-1, C2] U T) ^ C 2x2 such that 



(RH2.2) Y has continuous boundary values on C \ {— 1 , C2 } , denoted by Y + and Y_, such that 
Y + {z) = Y-(z)Jy(z), where 



Jy(z) 




zeT: 



z G 



-U2) 



with 



7 def -. 

oL = 1 



-2Aniri 



lie 



-Arnri 



sin(^4n7r) , 



(4.2) 



and with w{z) = w(z;A,B) defined in (|2.11j) . 
(RH2.3) As z -» 00, 



Y(z) 



I + O [ - 

z 







(RH2.4) Y is bounded in a neighborhood of z = — 1 and z = £2- 

In both cases C \ C has two connected components, one containing z = 1 and the other 
containing infinity; we denote these components by Oi and Qqo, respectively (see Figure |HJ). 

The steepest descent analysis, that we are going to carry out next, introduces new contours 
which are unions of a finite number of curves and arcs on C. In order to simplify notation 
we will call all the end points and points of self-intersection of such curves singular points, 
and the rest will be regular points of the contour. Hence, we could rephrase (RH1.4) and 
(RH2.4) saying that Y is bounded in a neighborhood of all singular points of C. 



4.2 First transformation Y 1— > U 

In order to shorten notation, we use the Pauli matrix 

,, f A 



0-3 







-1 
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and denote x°" 3 
function 



x 
x~ 



i I . Also for the sake of brevity, it is convenient to introduce the 



H(z) = G(z)w(z) , (4.3) 
analytic and single-valued in C \ (S U (— oo, 1]). Using H2.6[) - ()2.12l) . we immediately get that 
• For the case C.2, 



H(z) = exp 



A + B + 2 



t 2 -i 



dt 



for z e c \ (ru (-oo,l]) . 



(4.4) 



• For the case C.3, 



H(z) 



exp 



exp 



A + B + 2 



A + B + 2 



C2+*0 



t 2 - 1 

R(t) 



dt 



e -niA exp 



t 2 -l 



dt 



A + B + 2 



t 2 -i 



dt 



forz G C\(Pc(E) U (-oo,l]) ; 
forzeInt(Pc(S))nC+, 
forz G Int(Pc(S)) nC" . 



Observe that the same convention as in (|2.9|) for the path of integration applies: 

lim H(z) = 1 . 

C+\E32-<2 

Furthermore, taking into account ()2.26j) . in the case C.3, 



(4.5) 



lim H(z) = e 

-\Pc(£)9£-<2 



(4.6) 



In both cases C.2 and C.3, the sets of trajectories N and J\f r , (r > 0), introduced in subsection 
12.31 may be characterized by the conditions \H (z)\ = 1 and \H (z)\ = e r / 2 , respectively. Figure 
|7| shows also the domains where |-ff(^)| < 1. 

Function H has continuous boundary values at regular points of SU(— oo, 1], which satisfy: 



In case C.2: 



H+(z) 



■.niA 



H-(z), 



zeT, 
ze (-1,1) 
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IH(z)l<1 





Figure 7: Regions where |-H"(z)| < 1 for cases C.2 (left) and C.3. 



• In case C.3: 



H+(z) 



Hl\z), 

e- niA HZ 1 (z) 

e niA H_(z), 



z6l, 

(Ci,C2), 
ze (-1, l)\[Ci,Ca], 



(4.7) 



e -Ki(A+B) H_(z), z G (-oo, -1) . 
This allows us to express the boundary values of G at the regular points of E in terms of H: 



G+(z) 
G.(z) 



Hl(z), and G + {z)G~(z) 



w±(z) 



(4.8) 



Now we are ready to introduce the first transformation of the RHP, with the aim to 

1 /2 

normalize it at infinity. For d n in (|4.2j) . let us fix any value of d n , and define 



Hi na3 Y(z)G(z)- nas , 



in case C.2 ; 



U ^ \d~ a3/2 K na3 Y(z)G{z)- n(7s dn 3/2 , in case C.3. 



(4.9) 



with At given by (|2.1UI) . Obviously, matrix U solves now a new Riemann-Hilbert problem. 
Taking into account (|4.8j) we can state it as: 



(RH3.1) U is analytic on C \ C. 
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(RH3.2) U has continuous boundary values at the regular points of C, denoted by U+ and £/_, 
such that U + (z) = U-(z)Ju(z), where, for (A, B) in case C.2, 



Ju = < 



and for (A, B) in case C.3, 



'H- 2n {z) 1 

Hl n (z) 
1 H 2n (z)\ 
1 J ' 



z e r, 

z G 7 + U 7" 





: 2n (z) 









H 2n (z 




: 2 "(z) 


1 









1 




)• 





1 





z £ (Ci,C 2 ); 

(-l,Cl)- 



(RH3.3) As z oo, 



(RH3.4) Matrix U is bounded in a neighborhood of the singular points of C. 
4.3 Second transformation U i— > T 



(4.10) 



(4.11) 



By ()4.1Uj) - (|4.11|l . the jump matrix J[/ has oscillatory diagonal entries on E, along with 
exponentially decaying (as n — > oo) off-diagonal entries elsewhere and away from £12 (see 
Figure [7J). The aim of the next step is to transform the jump matrices with oscillatory 
diagonal entries into matrices asymptotically close to the identity matrix or to matrices with 
constant jumps. To this end, we take advantage of an appropriate factorization of Jjj and 
"open the lenses" around contours C. 

In case C.2, we use the following factorization of the jump matrix for z £ T (where we 
have taken into account that H + = 1/H- on V): 



H 



1 

-2n 





(z) 1 



1 
-1 



1 







Hf n {z) 1 



(4.12) 



Thus, the problem of the oscillatory diagonal entries of the jump matrix for z £ T may 
be solved by opening the lenses around T as it is shown in Figure El (left). The new contours 
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Figure 8: Contours C T for T. 



Yl and Tr are also oriented from Q\ to C2, an d this gives us two new bounded regions, £Ir 
and Ql, as well as modified domains Qj = Q\ \ VLl and f2^, = Q.qo \ Qr; we also denote 
C T = C U Ti U Fr, with the orientation shown in Figure EJ left. 

Hence, taking into account ()4.12j) . we define in case C.2 a new matrix-valued function 
T:C\C T — >C 2x2 by 



T(z) = U(z) 



1 

H- 2n (z) l) ' 

1 o N 

-H- 2n {z) 1 



U(z) 



for z G U 
for z G fix, , 

for 2: G iljj • 



T 

00 > 



(4.13) 



It solves the following RHP: 
(RH4.1) T is analytic for z G C \ C T ; 

(RH4.2) T(z) possesses continuous boundary values at regular points of C T , T+ and T_, related 
by the following jump conditions: 



T+(z) =T^(z)J T {z), zeC 1 



RIEMANN-HILBERT ANALYSIS FOR JACOBI POLYNOMIALS 



24 



where the jump matrix Jt is 



J T (z) 



1 H 2n (z)\ 

1 ) 

1 

H- 2n (z) 1 



z e r, 

z S 7 + U 7~ , 

z e r L u r R . 



(RH4.3) T(z) has the following behavior at infinity: 

T(z) = I + 0(l/z) asz^oo. 

(RH4.4) T(z) is bounded in a neighborhood of the singular points of C T . 

In principle, we could take advantage of factorization ()4.12p also in the case C.3. However, 
the geometry here is more complicated; this procedure would eventually yield a constant jump 
on whole X, which has now two components. In order to give a unified treatment to both 
cases in the next step, we use now a different factorization for Jjj: 







d 



-i 



-d n HZ 2n (z) 



1 

d n Hf n (z) 1 



1 







-2Anni 



HZ 2n (z) 1 



o r 
-i o. 



Hl 2n (z) 



z e r, 
ze (&,&)• 



These factorizations suggest to open lenses in the way shown in Figure |HJ right, which yields 
the new contour C T , splitting C into domains ttj , O^, Q^, and Qr, as shown. Now we define 
the matrix- valued function T : C \ C J 



T(z) = U{z) 



C 2x2 in the following way: 

z en R , 



1 

-d n H- 2n (z) l) ' 

d- 1 \ 

-d n H- 2n {z)j ' 

I 0\ 
-H~ 2n {z) I J ' 

1 
r 2Ami H- 2n (z) 1 



z 6 Of , 
zG 

z e , 



(4.14) 
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Then T(z) solves the problem (RH4.1)-(RH4.4), with Jjj replaced by 

(Ci,C 2 ) 



J T (z) 




4.4 Construction of the parametrices 

Now, we can see (cf. Figure [7J) that we have a single open arc joining the branch points (r in 
case C.2, and (C11C2) in case C.3) where the jump matrix Jt is constant, and at a positive 
distance from these arcs, Jt is asymptotically exponentially close to the identity matrix. 
Hence, by ignoring the "close-to-identity" jumps and condition (RH4.4) we are lead to the 
following problem: find an analytic matrix-valued function N(z) = I + 0(1 /z), z — > 00, and 
having the jump 

N + (z)=N^(z) 

on r (in case C.2) or on (Ci ; C2) m case C.3, with the orientation "from (j to £2" chosen. A 
solution of this model RHP, which is not unique in general, is (cf. [U Ch. 7]): 



N(z) 



a(z)+a(z) 1 a{z)—a(z) 1 \ 
2 2l I 

a(z)-a(z)- 1 a(z)+a(z)- 1 I ' 
2i 2 / 



(4.15) 



where a has been defined in ()2.14j) : it satisfies 

N(z) = 0(\z-Cj\- 1/4 ), j 



1,2. 



showing that the singularities at Q are L -integrable. Observe that the (1, 1) and (1, 2) entries 
of N coincide with Nn and N\2, introduced in (|2.1«-S|) . 
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We may expect N to be close to T away from (j and (2- However, in a neighborhood of the 
branch points the ignored jumps are no longer close to identity, and a different parametrix 
(model problem) is required. Now we look for two matrices P^\ j E {1)2}, which have 
the same jumps as T in a neighborhood of z = Q, and match TV" on the boundary of these 
neighborhoods. 

The construction of these matrices is well described for instance in [I]. Denote by A e (s) = 
{z € C : \z — s\ < e}, where we take e > sufficiently small. A local parametrix P^) in 
Ae(Cj)> j £ {1, 2}, solves the RHP with the same jumps as T there (see Figure EJ): 




Figure 9: Local analysis for cases C.2 (left) and C.3. 



(RH5.1) is analytic for z E A e (Cj) \ C T ■> bounded and continuous in A e ((j) \ C T ; 

(RH5.2) pO")(z) possesses continuous boundary values at regular points of C T Pi A £ ((j), P^ and 
(i) 

P_ , related by the following jump conditions: 

p (j) (z) = pV)( z )j pU) ( Z ), z£ c T n A £ (Q . 

(RH5.3) there exists a constant M > such that for every z G <9A £ ((j) \ C T , 

\\P^(z)N- 1 (z)-I\\ < — . 

n 

We describe the construction for £2; in order to simplify notation we write P instead of P^ 
whenever it cannot lead us into confusion. The jumps Jp = J P (2) specified in (RH5.2) are 
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In case C.2: 



Jp(z) 



In case C.3: 



Jp(z) 



1 
-1 Oj ' 
1 

H- 2n (z) 1 
1 H 2n (z)\ 
1 J 



o r 
-i o, 



i 



e- 2A ™F- 2n (z) 1 

1 o\ 
H- 2n (z) l) ' 
1 e 2A ™F+ 2n (z)\ 
.0 1 J 



zernA £ (c 2 ), 
j , zg (r L ur fl )nA £ (( 2 ) 
, z G 7" n A £ (C 2 ) • 



z G (Ca-e,<2), 

zG7 2 + nA £ (c 2 ) 
z g 7 2 " n A £ (c 2 ) 



z G (C 2 ,C2 + e) 

z g r. 



In order to solve the Riemann-Hilbert problems for P, let us first make a simple change 
of functions yielding piecewise constant jump matrices. For this purpose, we set for z G 
A^C2)\C T , 



R(z) = P(z) 



e -n<p(z)a^ in caS e C.2, 

e Anina 3 e -n<t>{z)o 3 j n cage q a 



(4.16) 



where <j) is the function introduced in (|2.2H) . In order to compute the new jumps we need to 
find how cp is related to H. In the case C.2, by (jOJ), exp(-</>(z)) = H{z) for z G A e (C 2 ) \ T. 
In the case C.3, by continuity of 4> in A e (£ 2 ) \ ((j, £ 2 ), 



exp(-c/>(z)) = < 



zg (fi L uo fl )nA e (( 2 )nc+ 



z g nj n A £ (c 2 )nc+ 



ff(z), 
^(z), 

e-^H-^z), zeO[nA £ (( 2 )nC-, 

e— iA tf(z), zg (fi L u^)nA £ (( 2 )nC- 



(4.17) 



and 



(0+ + (z) = 2niA, zG(C 2 -e,C 2 ). 
Now we can compute the jump matrix for R: J/j = e n <t>-( z ) (T 3 j T e -n<j} + (z)a 3 ^ name ly : 
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• in case C.2, 



Jr{z) 



= < 




z g r n A £ (c 2 ) , 



z e (r L u r fl ) n A £ (c 2 ) , 



* G 7- n A e (C 2 ) . 



• in case C.3, 



Jr(z) 



< 




) 



zG (Ca-e.Ca) 



z G 7 2 n A £ (C 2 ) 



2 G (C2,C2 + e)- 



Observe now that we have essentially the same local problem in both cases; in this way, we 
have reduced the RHP to the one studied in jH] (see also [U Ch. 7]), and we can write its 
solution explicitly 

For e > small enough, function 



defined in (|2.'22|) . is a conformal mapping from the neighborhood of the branch point onto a 
neighborhood of 0. In case C.2, T and 7~ are mapped onto the negative and positive real 
axis, respectively, and in case C.3 it happens to (Ci, C2) and (£ 2 , 1). Also we may deform the 
other curves (T^ and Tr in case C.2, and 7^ in case C.3) in such a way that the points on 
their image by / close to the branch point have the argument ±27r/3. 
Then the problem for R is solved by 



/(*) = \ (Hz)) 2/3 



R(z) = * (n 2 / 3 /(z)) 
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where \E' is built out of the Airy function Ai (see e.g. pQ) and its derivative Ai' as follows: 
Ai(i) Ai(u 2 t) 



Ai'(i) cj 2 Ai / (w 2 t) / 
Ai(i) Ai(a; 2 t) 



e e 



<T3 



e e 



1 

-1 1, 



Ai'(t) w 2 Ai'(w 2 *) 

Ai(t) -w 2 Ai(u;t)\ _« CT3 /l N 

Ai'(i) -Ai'(cut) J 6 \1 1, 

Ai(t) -cj 2 Ai(wt)\ _™ 

Ai'(t) -Ai'(wt) 



e e 



o-.'i 



< argt < 2vr/3; 
27r/3 < argt < 7r; 
— 7r < argt < — 27r/3; 
-2vr/3 < argi < 0, 



(4.18) 



and ijj = e 27r */ 3 . Finally, matrix P solving (RH5.1)-(RH5.3) is 

< , e n0( z ) CT3) in case C. 2, 



P(z) = E(z)R(z) 



e -A-Kinaz e n4>{z)a 3 ^ j n cage q 3? 



(4.19) 



where the analytic matrix function E is 

«i / 1 -l\ (nWfizfl^ 



4.5 Final transformation T *—>■ S 



a(z) 



(4.20) 



Now we may use matrix valued functions iV and for the final transformation. Recalling 
the definition of the contour C T , define the matrix- valued function S: 



s(z) « J r(^)iv^)- 1 , z G C \ (C T U A £ (Ci) U A e (C 2 )) 



(4.21) 



It is immediate to check that S 1 is analytic in C \ C s , where C s is the contour shown in 
Figure ITU1 Moreover, S : C \ C s — ► C 2x2 satisfies the following Riemann-Hilbert problem: 

(RH6.1) S is analytic in C \ C s . 

(RH6.2) S has continuous boundary values at regular points of C s , denoted by S+ and S_, such 
that S + (z) = S-(z)Js(z), where 



Js{z) 



{p^(z)N(z)-\ zedA E ((j), i = 1,2; 
{N( Z )J T (z)N(z)-\ zeC s \ (9A e (Ci) U <9A e (C 2 )). 
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Figure 10: Contours C s for S in cases C.2 (left) and C.3. 



(RH6.3) S(z) =I + 0(l/z), z^oo. 



(RH6A) S(z) is bounded in a neighborhood of the singular points of C s . 

Observe that by construction, J$ = I + 0{l/n) as n — > oo on dA £ ((j), j = 1,2, and is 
exponentially close to / on the rest of contours of C s . Using the same arguments as in [I] we 
conclude that 



uniformly for z £ C \ C . 

5 Proofs of the main results 

We establish strong asymptotics for {p n } tracing back all the previous transformations. For 
the sake of brevity, we do it explicitly only for the case C.3. The proofs in the case C.2 are 
very similar. 

5.1 Proof of Theorems [T] and El in the case C.3 

Assume that -1 < A < < B (case C.3); by 




Y{z) = d^ /2 K~ na3 U{z)G(z) na3 d- a3/2 



so that 




(5.1) 
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Assume that z £ C \ Pc(£), away from the branch points; without loss of generality we may 
take ze!l^ \ (A e (Ci) U A e (C 2 )) (see Figure IH1 or gO}) . Then by (|4l4~|) . 

U(z) = S(z)N(z) , 

and taking into account the expression of N in (|4.15ft . we obtain that uniformly on compact 
subsets of 0^ \ ( A^CiT U A7IC^) ) , 

p n (z) =Y 11 (z) = (^) n (SN) n (z) = (^) n N u (z) (l + O^j , (5.2) 

which proves ([2.15ft . 

If z belongs to the bounded component of C\T, we may assume without loss of generality 
that z e Uj \ A e (C 2 ). By KW\ and (jHHJ), 

U(z) = S(z)N(z) 
and again uniformly in compact subsets of Qj \ A e (^2 



H- 2n {z) -d- 1 
d n 



p n (z) = (^-T {[SN) n (z)H- 2n (z) + d n [SN] 12 (z)) 



G(Z)\ ( H -2n {z)Nii{z) ( 1 + (l\\ + dnNl2{z ) (l + of 1 



k J \ \ \n J J \ \n 

which proves ()2.16|) . Obviously, this formula is valid also if z £ T_, that is, if z lies on the 
"—"-side of r, away from the branch points. 

Assume now that z S T + away from £2. Again, without loss of generality we may take 
zeO K \ A e (C 2 ). By (jHH) and (jOT]) . 



and uniformly in compact subsets of £Ir \ A e (^2), 

Pn ( z ) = f^lj ([SN] n (z) + d n H- 2n {z) [SN] 12 (z)) 

G{Z) ^ Un(z) (l + O^jj + d n H~ 2n (z)N 12 (z) (i + oQ 
which proves ()2.18|) . Using (|4.7|) it is easy to see that formulas (|2.16|) and (|2.18|) match on V. 
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Finally, if z lies on the ±-side of the interval (£1, C2) 5 we assume z € fi^\(A £ ((i)UA £ ((2)), 
where by fUlj l and (jOJ) , 



U(z) = S(z)N(z) 




z G Q L 



Hence, uniformly in compact subsets of fij \ (A e (d) UA £ (^)), 
Pn(z) = (^) n ([SN} U (z) + H-^{z) [SN] 12 (z)) 

l + 0^+H-^(z)N l2 (z) (l + oQ 



while uniformly in compact subsets of Q L \ (A E (£i) U A £ ((^2)), 
Pn(z) = (j^-J ([•SWlii (2) " c- 2A '"'H- 2 "(z) [SN\ 12 (*)) 

;^m)" ( WllW (i +0 (i))_ e -«*ir* WWuW (i +0 (i))) 



This finishes the proof of (t2~TT)l) - (l2~2fil) . 

5.2 Proof of Theorem |U 

Let z E A e (C 2 ); then by (JOT}, 

[/(z) = 5(z)P(z)^- 1 (z), 

where if (z) is one of the matrices given in the right hand side of (|4,14|) . Gathering (|4.19|1 - 
(jOOJ), we get that 

with i„ = f n 2 / 3 f(z) and * given by (|4.18|) . For instance, if z € A e (£ 2 ) n Or, using (|4.17|) we 
get 

. . / 1/4 \ °"3 



x g— Amnas e n4>(z)a3 



a(z) ^ 
1 ' 0\ 
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Observe that for z G A E (£2) n 0^, we have < aigf(z) < 2ir/3, and we use the expression 

*(*) 



Hence, 



(Ai(t) Ai(cu 2 t) 
\Ai'(t) uj 2 Ai'iojH) 



where 

Consequently, 



U2i{z) ) \-% -%) \a(z) / \A{t r 



A(t) = e- Anin Ai(t) + 2» eT sin(^n) Ai(u 2 t) , w = e 2 ™ /3 



1 + Of - 

n 



trilW = ^»(^^) (i +0 (i))-^^ 



and taking into account ()5.1j) and the fact that in A e (^2) H fi^, exp(— 0) = ff, we arrive at 
Q2.23JI . Proceeding in a similar way, we see that this expression is also valid for z in the other 
regions of A £ (^2)- 

5.3 Proof of Theorem [5] 

This theorem is a corollary of Theorems ^ and [21 First, taking into account (|2.15|) and that 
function Nn, defined in ()2.13|) . has no zeros in the plane cut from (j to C2, we see that zeros 
of {p n } cannot accumulate at C \ Pc(S). 

Consider in particular case C.3. Now the asymptotic location of the zeros of Jacobi 
polynomials depends also on the value 

e~ r = lim \d n \ 1/n = lim | smOrn)| 1/n , 

(assuming it exists), where d n , defined in (|4,2|) . depends upon the distance of a n = An to 
the integers, in such a way that 

e~ r = lim (dist(a n ,Z)) 1 / n . 

n— >oo 

Let z S Int(Pc(S)), that is, z lies in the bounded component limited by the contour V. We 
can choose e > small enough such that z ^ A e (£ 2 )- Then, the asymptotic formula l|2.16j) . 



Pn (z) =-L ((G(z)w 2 (z)) n N n (z) (l + Q 

+2ie~ An7Ti sin(Avrn) G n {z)N 12 {z) (l + o(^ 
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is valid. Hence, z is a zero of p n only if 



H- 2n (z) = -2ie- An7Ti sinMvm) (l + o(- 

N u (z) V \ n 

Since N^/Nxi is uniformly bounded and uniformly bounded away from zero, we see that the 
zeros in this domain must satisfy 



Itf- 1 ^)! = | sinOm)! 1 /^) (l + O Q jj 



e -r/2 (i + (I 
n 



It remains to use that = e r / 2 defines in Int(Pc(E)) the curve r r . 

Once we have established where the zeros accumulate, it remains to prove that they 
asymptotically distribute according to the corresponding measures in parts (i)-(ii). To this 
end we can use the second order linear differential equation satisfied by Jacobi polynomials 
y n = Pt M (see e.g. (33 §4.22]): 

(1 - z 2 ) y"(z) + [f5 n -a n - (a n + P n + 2)z] y' n (z) + n(n + a n + n + 1) y n (z) = . 

If we rewrite it in terms of h n = y' n /(ny n ), we obtain a Riccati differential equation: 

M 2x ( 1 , / / s . ,2 / x\ , Pn ~ a n ~ (a n + Pn + 2)z , U + a n + /3 n + 1 , . 

(1-z -h n {z) + h n (z) + K\ z ) + = °- ( 5 - 3 ) 

\n J n n 

Let v n denote the normalized zero counting measures of y n = P^ n ^ n \ By a weak compact- 
ness argument we know that there exists an infinite subsequence A C N and a unit measure 
v such that v n — > u, n £ A, in the weak*-topology. In the first part of this proof, we saw 
that supp(^) consists of a finite union of analytic arcs or curves, and every compact subset 
of C \ supp(i/) contains no zeros of pjf"'^ for n sufficiently large. 
Hence, 

dv n {t) f dv{t) 



h n {z) = \ — ^--r — ► H z ) = / — ~7 1 »eA, 
J T z — t Jy z — t 

locally uniformly in C\supp(^). Taking limits in (|5.3|) we obtain that h satisfies the quadratic 
equation 

(1 - z 2 ) h 2 {z) + [B - A - (A + B)z] h{z) + A + B + 1 = , 

so that 

f dvit) A + B + 2 R(z) 1 / A B \ 

/ —4 = o 2H " 9 7 + — T > z G C \ supp " • 

jYZ — t 2 z z — 1 2 \z — 1 z+lj 

By Sokhotsky-Plemelj's formulas, on every arc of supp(^), 



, , , A + B + 2 R + {z) 

MZ) = 2 m z 2 -l dZ - ^ 
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Now, we are concerned with proving part (ii) of the theorem, related to the case C.3. First, 
consider the generic case when r = 0. In this case, the measure \x in Q2.5JI is supported on 
£ = r U [Ci, £2]- Thus, by l|5.4|) . //' = u' a.e. on supp(z/), /u, v being unit measures. Therefore, 
v = fjL. The proof in the case < r < 00 is similar, but with measures fj, r , given in l|2.24|) . in 
place of fx. Finally, for the degenerate case r = 00, which takes place when parameters a n 
approach the integers faster than exponentially, it is enough to take into account that the 
Cauchy transform of the measure da = — A5\ is a(z) = —A/(z — 1). 

Finally, the proof for the Case C.2 is totally analogous. 
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